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Abstract 

The bilayer Heisenberg antiferromagnet is known to exhibit a quantum- 
critical transition at a particular value of the inter-layer coupling. Using a 
new type of coherent state, appropriate to the special order parameter struc- 
ture of the bilayer, we map the problem onto the quantum non-linear sigma 
model. It is found that the bare coupling constant diverges at the classical 
transition of Chubukov and Morr, so that in any finite dimension the actual 
transition occurs inside the ordered phase of the classical theory. 
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The study of non-classical collective quantum states of matter is a central theme of 
modern condensed matter physics. Despite the successes in 1+1 dimensions, it has proven 
difficult to address these matters in higher dimensions. Either the minus-sign problem 
intervenes (as in, e.g., the t-J model and frustrated spin models), or the tendency towards 
classical order is too strong (e.g., unfrustrated spin models). The class of bilayer Heisenberg 
models is special in this regard HIQ]. It is sign-free, and convincing numerical evidence 
exists showing that its long- wavelength behavior is governed by the 0(3) quantum non- 
linear sigma model (QNLS) with tunable bare coupling constant u ||. The relationship 
between the microscopic model and its long wavelength behavior is non-trivial. Chubukov 
and Morr (CM) made the key observation that, in order to construct the classical limit, the 
severe local (interplanar) fluctuations have to be integrated out first ||. In the resulting 
singlet-triplet representation, a phase transition between a Neel state and an incompressible 
state is found already at the classical level. CM conjectured that this transition corresponds 
to the quantum critical transition found in numerical studies. Here it is shown that this 
is not correct. Because of the special structure of the order parameter, the standard su{2) 
generalized spin coherent state does not suffice for the construction of the path integral. We 
introduce a novel type of coherent state which allows us to straightforwardly recover the 
QNLS describing the long wavelength behavior. We find that the bare coupling constant 
of the field theory diverges at the classical transition of CM. The quantum phase transition 
therefore occurs well before the classical transition can occur, and the latter is therefore in 
any finite dimension an artefact. 

It is convenient to consider the "bilayer" model in arbitrary dimensions, with an added 
magnetic field (B), 

TC — J± ^ (Sjl ■ Sji + S i2 ■ Sj2) + 

<ij> 

+Ji^2s a • s i2 - B + s; 2 ), (1) 

i i 

where < ij > runs over the bonds of two <i-dimensional hypercubes 1 and 2. The anti- 
ferromagnetically coupled {J\ > 0) s = \ Heisenberg spins s ir] are coupled locally by J 2 . 
Following CM, we first integrate out the J 2 term ||||. Define the sum and the difference of 
the spin operators, 

S = s i + s 2 , 

S = si - s 2 , (2) 

such that 

w = l J i E $ ■ 4- + ^ ■ + 

1 <ij> 

+^E(5?-^)-£-E& (3) 

Eq. (g) amounts to a transformation to a singlet-triplet basis. Introducing hard-core bosons 
creating the local singlet state, a\ = -^(c]ij. c i2t ~ c !iT c i2j,)' anc ^ ^ ne l° ca l triplet i,o,— l 
= c liT c i2T etc)) Eq. (0) can be alternatively written as, 
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S z = blh - 6Li6_i, 
S + = V2 (b\b + 6j6_! 
= -atfco - 6j 



a, 



5+ = \/2 (61a - a+fe-i) . (4) 

— * 

5 describes S 1 = 1 spins, while 5 is related to fluctuations from triplets to singlets. These 



(5) 
(6) 
(7) 

At the J 2 = point (two decoupled layers) the problem has an 0(4) global invariance, which 
is broken for any finite J 2 , leaving only an invariance under the SU(2) subgroup Eq. (^j). 
The unconventional aspect of this problem is that for positive J 2 the spontaneous symmetry 
breaking involves the generators S. The J 2 > classical saddle-point of CM is easily seen 
to correspond to the vacuum amplitudes (z = 2d), 



[S a , 


S b ] 


= V^abcS , 


[S a , 


S b ] 


= lEabcS , 


[S a , 


S b ] 


= l^abcS ■ 



n = i<E(- 1 ) <5 «> 



j. 



1 " -1T^ Q ( J i z - J 2) n, (8) 



"a = ^<E 4^) = ^ (1 + 5^) © ( - J*) + 

+e (j 2 - j lZ ) , (9) 

where h is a vector on the unit sphere. The condensation of S (Eq. ||) and the existence 
of a mean singlet density (Eq. ^) is a direct ramification of the explicit symmetry breaking 

due to the interplanar coupling. Q is still a vector order parameter, because S transforms as 
a vector under S. It is therefore a Neel order state, albeit one with a variable local moment 
size, which implies that its long wavelength behavior should be described by the QNLS. 

On the classical level it is found that is nonzero for all positive J2, while Q vanishes 
continuously at J 2 = J\Z, where becomes equal to one. This is the transition found by 
CM. Regarding its formal status, it is easily seen that this classical theory becomes exact in 
infinite dimensions |5|,0]. The energy fluctuations disappear in this limit: AE/E oc 1/yNd. 

In addition, we note that fl also exists in 2+1 D, at least in the vicinity of the quantum 
critical point: the correlation functions in terms of s\ and s* 2 go to zero at the transition 
with their ratios fixed according to Eq. (|j) 0. 

What is wrong with the assertion that this transition and the quantum-critical transition 
in 2+1 D are the same? The transition in infinite dimensions is a classical transition. In 
terms of the singlet-triplet basis, the quantum fluctuations disappear at the lattice cut- 
off and thermal fluctuations dominate at any finite temperature. The numerical study 
shows quantum criticality ||: at zero temperature, the quantum fluctuations are scale 
independent. In the remainder we will show that this classical theory becomes pathological 
in the neighborhood of the classical transition. 
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Coherent state path integrals offer a convenient framework to study quantum order 
parameter fluctuations ||. Because of the special status of the order parameter Eq. 
the usual generalized spin coherent states do not suffice. Our key result is the discovery 
of a special coherent state for this type of order parameter structure. Next to the general 
requirements of normalizability and the existence of the identity, it should be demanded from 
coherent states that they reproduce all properties of the classical sector. Besides reproducing 
Eq.'s they should also allow for an S derived vacuum expectation value, 

fi = -^<E3>- ( 10 ) 

We find that the following coherent state satisfies all these requirements. 

\QQ)=e^ sz e iesy e^ sx e^ v \ X ), (11) 

with the reference state, 

Ix) = (cosxa 1 " + sinx&o)|vac). (12) 
Eq. ( pi]) looks conventional. It refers to the various rotations related to the 0(4) 



symmetry. The novelty is Eq. (|i~2"D: instead of the usual maximum weight state, this 

non-exact state underlies the order parameter structure Eq.'s (|8],|9|), with Q chosen along the 
z axis, while x is fixed by the explicit symmetry breaking interaction ~ J 2 . The freedom 
implied by Eq. ( |TT| ) might at first instance appear as redundant. However, it turns out 

that the stiffness in the temporal direction is caused entirely by the fluctuations from Q into 
the Q direction, and the four angles appearing in Eq. ( |T2"D take care of the independent 

rotations of Q and Q. Explicitely, ip parametrizes a rotation from Q to Q _L Q (S ■ S = 0). 

The rotation of Cl in the plane perpendicular to Q is parametrized by 9 2 ■ This is the only 

free rotation left to Q in a magnetic field. 9 and fix the direction of Q. 

We obtain the following expressions for the vacuum amplitudes with respect to this 
coherent state 

riA = cos 2 x cos 2 (13) 
Q = sin 2x sin V(— cos 9 cos <f), — cos 9 sin </>, sin 9), (14) 

= sin 2% cos cos 9 2 9(9, 0) — sin # 2 0(0) , (15) 

where 9 and are the local unit vectors in the 9 and direction, 9 = 
(sin 9 cos 0, sin 9 sin 0, cos 9) and = (sin 0, — cos 0, 0). The identity becomes, 

i = Jdfi(n,n\ (16) 

9 rn/2 fir/2 r2n 

— d X d^ / d9 2 x 



7T 4 ^0 " JO JO 



x /^d0 [ W/2 d9cos9\Qtt)(tttt\ 

JO J-n/2 



4 



By taking expectation values with regard to \QQ) (classical limit), we find the 0(3) 
invariant version of the mean-field theory of Chubukov and Morr. Minimization of the 
classical energy with regard to the coherent state angles yields, 

cos2 X0 = ^ + O(£ 2 ), (17) 

^i/irJ^- (18) 

— * 

with 9 and (f> fixed such that Q points in the direction of the magnetic field. We recover the 
classical order- disorder transition at J 2 = JiZ, where both Cl and the induced magnetization 
Q vanish according to Eq.'s (|T^JT5|) . 

The derivation of the path integral is standard ||. Using the Trotter formula, the 
evolution operator in imaginary time is written as (Nt is the number of time slices, St the 
imaginary time interval, N t St = (3), 

Z = lim (e-**)"*. (19) 

it -> 

Inserting the identity ([T^) at every intermediary time and expanding Z to lowest order in 

St, 

N t 



Z= lim fVfif[ [(nf2(t,)|fifi(ti+i)> 

Nt -» oo J L 
A, — » n t—X 



-Stinn^imnn^))] , (20) 

where the integration measure T>fi is given by d/i ff2j,f2^, while {ti} is the set of 

intermediary times in the imaginary time interval [0,/3]. The kinetic term in the action 
follows from the first term inside the square brackets, 

(nn(t,)|nnfe + . 1 )) = (1 + o(s 2 )) e **fe>+ <«?>, (21) 

with 

$ = ^sin2xi sva , 4>i(sm9idt4>i + d t 9 2 i) 

i 

= -]Tsin2 X A-^SxS, (22) 

— * —* 

where = 0,/ sin2x and = 0,/ sin 2%, so O 2 + O 2 = 1. 
The potential energy is (B = 0), 

V = y £ sin 2 Xi sin 2 Xj (<5< ■ O j + { ■ + 

+ iE(l-40>s 2 x,). (23) 
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Taking the time continuum limit, the path-integral becomes Z = lim N t -> <x> / T>fxe Sm , 

«t -» 

with the real-time action, 

S M = [ T dx [-$(x ) + V(x )} (24) 
Jo 

To derive the long wavelength theory, O and O are separated into a slowly varying order- 
parameter part and a rapidly fluctuating part which will be integrated out. The fluctuations 
in x are massive because of the explicit symmetry breaking, and can be neglected. We are 
left with, 



Oi = T]i [mi + aL\\ij + aL ±i} (25) 
Oi = rhi + aLi. (26) 

The (staggered) fluctuation L\\i is parallel to the order parameter rhi (Vi = ±1 depending 
on the sublattice). L has a component along m, but is perpendicular to m because of the 

constraint Oi - Oi = 0. As we already indicated, despite the fact that the order-parameter 
part of O is zero in the absence of a magnetic field, the fluctuations in this quantity are 
actually producing the stiffness in the time direction and should be carefully integrated out. 
We expand to second order in the lattice constant a, which will be taken to zero at the end of 

the calculation. Using the constraint Of + Of = 1, the fluctuation Liu is eliminated from the 

action. Different from the single-layer system, two canting-fields result, L±i and Li, which 
have to be integrated out. The former does not influence the long wavelength behavior, 
while the latter, which is related to the response of the system to a uniform magnetic field 
applied perpendicular to the plane of ordering, is responsible for the stiffness in the time 
direction. 

After expanding in a and eliminating L\\, the kinetic term becomes 



$ = — sin 2\i-r^Li ■ doirii x rhi + stagg- terms. (27) 



Using rhi ~ m j — a di- t jrhi, it can be seen that the staggered terms give contributions which 
are of third order in a. The expression for $ is identical to that for the single-layer system, 
apart from the factor sin 2% and the absence of a topological term. Within the limitations 
of the semiclassical expansion, the above derivation is in principle valid for any dimension, 
including the 1+1 dimensional two leg spin ladder systems. The usual argument for the 
irrelevance of topological terms in these systems are based on the proximity of Neel order 
on both chains separately: the topological terms in the two rows cancel each other. Here 
we find that this holds regardless the strength of the local fluctuations. We notice that, 
according to Haldane's conjecture 0, the spectrum of the two leg ladder has to be gapped 
for any J 2 ^ 0. 

The potential term is written in the form 



V = Ji sin2xisin2xj 

<i,j> 



°- (d i ^ j m t ) 2 + 
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+a 2 (pi + L 2 ±l ) - 2] + J 2 J2 a 2 L 2 cos 2 X 

i 

J; 



4cos 2 Xi . (28) 



In the continuum limit (a — > 0), the summations over sites are replaced by integrations 
over space, Yh ~^ a ~ d I d d x. The 0(1) term in Eq. (p8"D, corresponding with the mean field 
energy for the bilayer model, acquires a large prefactor a~ d and can be integrated by steepest 
descent. This yields the mean-field expression for x, Eq. ( jlT]) . 

After integrating over the fluctuations L and Lj_ we recover the effective action respon- 
sible for the long wavelength fluctuations, which is the 0(3) QNLS 



Sm = \j d d+1 x X ± (domf -p s J2 (d, 



jm 

a=l 



(29) 




Although the form of Eq. (P5|) is dictated by symmetry, the parameters appearing in the 
effective theory have a quite different meaning in terms of the microscopic model than is the 
case in single layer problems. Taking the saddle point values, the perpendicular susceptibility 
and the spin stiffness become respectively, 

(30) 
(31) 

Both the susceptibility and the spin stiffness vanish at the classical transition at J 2 = zJ\. 
The spin-wave velocity v s = \Jp s /x± remains finite at the transition, and no divergencies 
occur on the Gaussian level [|]]. The stability of the classical state against quantum melting 
is, however, controlled by the dimensionless coupling constant u = a 1 ~ d / ^p s X± which is 
found to diverge at the classical transition as u ~ 1/^1; "where J 2 is the reduced interlayer 
coupling = (J\Z — J 2 )/^2- In any finite dimension, the 0(3) QNLS quantum critical 
transition occurs at a finite value of the coupling constant and it follows that the long wave- 
length fluctuations destroy the Q type Neel order before the classical critical point is reached. 
Accordingly, the quantum critical transition of the bilayer model is of the 0(3) QNLS kind, 
and the classical transition of Chubukov and Morr only exists in infinite dimensions. 

This theory is even quantitatively reasonable. One loop renormalization theory for the 
QNLS in 2+1 dimensions puts the critical coupling at u* = 4tt ||. Using the saddle point 
values for the spin stiffness and susceptibility (Eq. [H|), we find the quantum transition to 
occur at J%/Ji = 3.3. Given that l/S'-like corrections are neglected 0], the agreement with 
the value of 2.5-2.6 obtained from Quantum Monte Carlo and series expansions is 
reasonable. 

In summary, we have clarified the origin of the quantum critical transition of the bilayer 
Heisenberg problem. The key aspect is that the order parameter structure as discovered 
by Chubukov and Morr is unusual. Although this order parameter is macroscopically of 
the usual 0(3) vector kind, and therefore described by the 0(3) quantum non- linear sigma 
model, its microscopic status is unconventional. The operators acquiring a vacuum am- 
plitude (S) are not the ones expressing the global SU(2) invariance of the problem. This 
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kind of order parameter structure arises naturally in the present context and we expect it 
to be quite common in the general context of quantum magnetism [10|. Our main result 
is the discovery of a new type of spin coherent state which allows for the requantization 
of such order parameter structures. As applied to the bilayer problem, the novelty is that 
in any finite dimension the classical theory becomes highly pathological: the bare coupling 
constant of the field theory diverges at the classical transition, explaining why the quantum 
transition obeys 0(3) QNLS universality. 
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